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ABSTRACT 



Aims. We present the first multi-epoch study that includes concurrent mid-infrared and radio interferometry of an oxygen-rich Mira star. 
Methods. We obtained mid-infrared interferometry of S Ori with VUTI/MIDI at four epochs in December 2004, February/March 2005, 
November 2005, and December 2005. We concurrently observed v = 1, 7 = 1 - (43.1 GHz), and v = 2, / = 1-0 (42.8 GHz) SiO maser 
emission toward S Ori with the VLBA in January, February, and November 2005. The MIDI data are analyzed using self-excited dynamic 
model atmospheres including molecular layers, complemented by a radiative transfer model of the circumstellar dust shell. The VLBA data are 
' reduced to the spatial structure and kinematics of the maser spots. 
Results. The modeling of our MIDI data results in phase-dependent continuum photospheric angular diameters of 9.0 ± 0.3 mas (phase 0.42), 
7.9 ±0.1 mas (0.55), 9.7 ±0.1 mas (1.16), and 9.5 ± 0.4 mas (1.27). The dust shell can best be modeled with AI2O3 grains using phase-dependent 
inner boundary radii between 1.8 and 2.4 photospheric radii. The dust shell appears to be more compact with greater optical depth near visual 
minimum (tv ~ 2.5), and more extended with lower optical depth after visual maximum (tv ~ 1.5). The ratios of the 43.1 GHz/42.8 GHz SiO 
maser ring radii to the photospheric radii are 2.2 ±0.3/2. 1 ±0.2 (phase 0.44), 2.4 ±0.3/2.3 ±0.4 (0.55), and 2.1 ± 0.3/1.9 ± 0.2 (1.15). The maser 
spots mark the region of the molecular atmospheric layers just beyond the steepest decrease in the mid-infrared model intensity profile. Their 
^ velocity structure indicates a radial gas expansion. 
k>( Conclusions. S Ori shows significant phase-dependences of photospheric radii and dust shell parameters. AI2O3 dust grains and SiO maser 
spots form at relatively small radii of ~ 1.8 - 2.4 photospheric radii. Our results suggest increased mass loss and dust formation close to the 
d surface near the minimum visual phase, when AI2O3 dust grains are co-located with the molecular gas and the SiO maser shells, and a more 
expanded dust shell after visual maximum. Silicon does not appear to be bound in dust, as our data show no sign of silicate grains. 

Key words. Techniques: interferometric - Masers - Stars: AGB and post-AGB - Stars: atmospheres - Stars: mass-loss - Stars: individual: 
SOri 
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* Based on observations made with the Very Large Telescope 
Interferometer (VLTI) at the Paranal Observatory under program IDs 
074.D-0075 and 076.D-0338. 

** Based on observations made with the Very Long Baseline 
Array (VLBA) under project BB192. The VLBA is operated by 
the National Radio Astronomy Observatory (NRAO). The National 
Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associate 
Universities, Inc. 

*** Tables ??-?? are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
^http ://cds web. u- s trasbg. f r/ abstract.html| 

^ Color versions of Figs. 2-10,11,12 are available in electronic form 
via |http://www.edpsciences.org| 



1. Introduction 

The evolution of cool luminous stars, including Mira variables, 
is accompanied by significant mass loss to the circumstellar en- 
vironment (CSE) with mass-loss rates of up to lO"'* Mo/year 
(e.g.. Jura & Kleinmann 1 1 9901 1 . This mass-loss process signif- 
icantly affects any further stellar evolution and is one of the 
most important sources for chemical enrichment of the inter- 
stellar medium. The detailed nature of the mass-loss process 
from evolved stars, and especially its connection with the pul- 
sation mechanism in the case of Mira variable stars, is a matter 
of cuiTent investigation. 

The conditions close to the stellar surface can be studied 
well by means of optical long-baseline interferometry. This 
technique has provided information regarding the stellar pho- 
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tospheric diameter, effective temperature, center-to-limb inten- 
sity variation (CLV), and atmospheric molecular layers for a 
number of Mira variables (see, e.g., Quirrenbach et al. 1992; 
Haniff et al. [T995l van Belle et al. |T996l Thompson et al. 2002; 
Mennesson et al. 120021 Hofmann et al. 120021 van Belle et al. 
12002 ^, Perrin et al. l2004l MiUan-Gabet et al.l2005T). 

Woodruff et al. (|2004] | and Fedele et al. (2005 ) have re- 
cently shown that observed near-infrared /iT-band visibilities 
of the oxygen-rich prototype Mira variables o Cet and R Leo 
are very different from uniform disc (UD) models already 
in the first lobe of the visibility function, and that they cor- 
respond closely to predictions by self-excited dynamic Mira 
model atmospheres that include effects from molecular lay- 
ers (models by Hofmann et al. [19981 Tej et al. l2003bl Ireland 
et al. I2004ai l2004b ). Ohnaka et al. ('2006b') studied the com- 
parison of the same dynamic Mira model atmospheres with 
mid-infrared interferometric and spectroscopic observations. 
Recently, Ireland & Scholz (2006) added the formation of dust 
in a self-consistent way to the same dynamic Mira model atmo- 
spheres. They find that dust would form at 2-3 times the aver- 
age photospheric radius for certain plausible parameter values. 

The structure of the atmospheric molecular shells located 
above the photosphere, as well as the dust shell, can be 
probed by mid-infrared interferometry. This has been success- 
fully demonstrated by Ohnaka et al. (2005) using the spectro- 
interferometric capabilities of the VLTl/MlDl facility for ob- 
servations of the oxygen-rich Mira star RR Sco. The model 
used in this work includes a warm molecular layer consist- 
ing of SiO and H2O, as well as an optically-thin dust shell 
of AI2O3 and silicate. Recent VLTl/MlDl observations of the 
carbon-rich Mira star V Oph by Ohnaka et al. (I2007I I indi- 
cate that carbon-rich Miras also have extended atmospheric 
layers of polyatomic molecules (C2H2) and dust shells (amor- 
phous carbon and SiC). Information on dust shells around Mira 
variables has also been derived using mid-infrared interfer- 
ometry with the Berkeley Infrared Spatial Interferometer ISl 
(see e.g., Danchi et al. 1994; Weiner et al. 120061 Tatebe et al. 
I2OO6I 1. Ireland et al. (I2005I I recently used Q.9 fim interferomet- 
ric polarimetry measurements of the Mira variables R Car and 
RR Sco to place constraints on the distribution of the dust shell. 

Complementary information regarding the molecular shells 
around oxygen-rich AGB stars can be obtained by observing 
the maser radiation that some of these molecules emit. Maser 
emission from the three most common maser molecules, SiO, 
H2O, and OH, traces regions of the CSE on scales from a 
few to several hundred AU. Masers provide a unique probe of 
structure, kinematics, and polarization properties of the envi- 
ronment of these stars. The SiO maser radiation in the CSE 
of oxygen-rich Mira variables has been mapped using the 
Very Long Baseline Array (VLBA) by, e.g.. Diamond et al. 
(fT994] |. Kemball & Diamond (fT997i i. Boboltz et al. (,1997 ). 
Diamond & Kemball (l2003l l. Cotton et al. (|2004] |. Boboltz 
& Wittkowski (l2005l hereafter BW05), Boboltz & Diamond 
(1200 5 ). and Cotton et al. (12006V The SiO maser emission is 
typically found to arise from a clumpy ring within a few stellar 
radii of the photosphere, indicating a tangential amplification 
process. Humphreys et al. ( 1996ll2002l l have predicted the lo- 
cation at which SiO molecules form in the CSE of Mira vari- 



ables by combining maser models with hydrodynamic pulsa- 
tion models. Their model predictions are roughly consistent 
with the measurements mentioned above. 

Observational results regarding the detailed relationships 
between the stellar photosphere, the molecular layer, the dust 
shell, and the SiO maser ring often suffer from uncertainties in- 
herent in comparing observations of variable stars widely sepa- 
rated in time and stellar phase, as discussed in BW05. To over- 
come this limitation, we have established a program of con- 
current infrared interferometry using the VLTl and radio inter- 
ferometry at the VLBA. The former aims at constraining the 
photospheric radius, the characteristics of atmospheric molec- 
ular layers, and parameters of the dust shell. The latter aims at 
concurrently mapping the SiO maser emission. Our final goal is 
a better understanding of the mass-loss process and its connec- 
tion to stellar pulsation. Our pilot study in BW05 on the Mira 
variable S Ori included coordinated near-infrared ^T-band in- 
terferometry to constrain the stellar photospheric diameter and 
VLBA mapping of the SiO maser radiation toward this source. 

A further uncertainty in comparing photospheric radii to the 
extensions of the dust and maser shells arises from the compli- 
cation that near- and mid-infrared CLVs of finite bandwidth 
include a blend of intensities from continuum-forming layers 
and overlying molecular layers. This effect has often resulted in 
over-estimated continuum photospheric diameter values (see, 
e.g. the discussions in Jacob & Scholz 120021 Mennesson et al. 
I2OO2I van Belle et al. l2002l Woodruff et al.'2004', Ireland et al. 
I2004ail2004 b; Perrin et al. 12004. Fedele et al. ,2005 j. An over- 
estimated photospheric diameter would result in biased relative 
distances of the dust shell and the maser ring from the stellar 
surface, even if obtained at the same stellar phase and cycle. A 
detailed comparison of observations to dynamic model atmo- 
spheres, as mentioned above, can be used to relate the observ- 
able quantities to the continuum photospheric radius and thus 
to overcome this limitation. 

Here, we present VLTI/MIDI mid-infrared interferometry 
of S Ori at 4 epochs/stellar phases and coordinated VLBA map- 
ping of the SiO 42.8 GHz and 43.1 GHz maser transitions at 3 
epochs/stellar phases that are contemporaneous to the first 3 of 
our 4 MIDI epochs. The first two epochs are located near the 
stellar minimum, and the later epochs shortly after the follow- 
ing stellar maximum. 

2. Lightcurve and characteristics of S Ori 

S Ori is a Mira variable star with spectral type M6.5e-M9.5e 
and V magnitude 7.2-14.0 (Samus et al. 2004). Templeton et al. 
(120051 ) report that S Ori's period appears to vary in a seemingly 
sinusoidal fashion between about 400 and 450 days over about 
the past 100 years. The mean period is P = 420.0 days with a 
latest period, corresponding to our observation epochs, of P w 
430 days. We find a good correspondence with the AAVSO 
(Henden et al. |2006] ) and AFOEV (CDS) data for the cycle of 
our observations, as well as for the most recent ~ 5 cycles using 
P - 430 days and a Julian Date of last maximum brightness 
T() - 2453190 days. Figure [T] shows the visual lightcurve of 
S Ori as a function of Julian Date and stellar phase based on 



M. Wittkowski et al.: The Mira variable S Ori: Photosphere, molecular layer, dust shell, and SiO maser shell 



3 



Stellar cycle/phase 
-3 -2 -1 



10 - 

- ^ 

12 - 

- ^ 

14- 
1000 




1500 



0.0 



2000 2500 3000 
Julian Date - 2450000 

Stellar ohase 
0.5 



4000 



8 

10 - 



1.0 


1.5 




AAVSO + 








AFOEV X 




+ 







^ 12 - 
> 

14 -- 



VLBA observations t t 
1g_ MIDI observations t 1t 
Epoch A B 



t 

t t 

C D 



3000 3200 3400 3600 3800 
Julian Date - 2450000 



4000 



Fig. 1. Visual lightcurve of S Ori as a function of Julian Date 
and stellar cycle/phase. Data are from the AAVSO (Henden 
et al. l2006l l and AFOEV (CDS) databases. We adopt a pe- 
riod of f = 430 days and JD of last maximum brightness 
To - 2453190. The arrows indicate the dates of our VLTI 
and VLBA observations, a) The most recent 5-6 cycles illus- 
trating the correspondence of the data with our adopted P and 
Tq values. The indicated dates of observations are the epochs 
from BW05 (JD between 2600 and 2800) and the epochs of 
the present paper (JD between 3300 and 3800). b) Enlarged 
view covering our VLTI/MIDI and VLBA observations pre- 
sented in this paper. The dates of observations are combined 
into 4 epochs A, B, C, D. 



these values. Also indicated are the dates/epochs of our VLTI 
and VLBA observations as described in Sects. [3]and|4] 

The distance to S Ori, as for many Mira stars, is not well 
known. Van Belle et al. (I2002I I have estimated the distance to 
SOrito480pc + 120 pc based on a calibration of the period- 
luminosity relationship by Feast et al. (il98^ . and we adopt this 
value. 

The near-infrared /T-band UD angular diameter 0^^ of 
S Ori has been measured by van Belle et al. ( 1996), Millan- 
Gabet et al. (.2005) . and by BW05 to values between 9.6 mas 
and 10.5 mas at different phases. Diameter measurements of 
S Ori at other wavelengths have so far not been reported. 

S Ori exhibits SiO and OH maser emission (Benson et al. 
119901 1. while a detection of H2O maser toward S Ori has not 



been reported. The maps of the SiO maser emission in BW05 
are the first toward S Ori. They derived average distances and 
standard deviations of the SiO maser spots from the center of 
their distribution at phase 0.73 of 9.4+ 1.4 mas and 8.8 + 1.7 mas 
for the 43.1 GHz and 42.8 GHz transitions, respectively. 

Sloan & Price (1998) report on a relatively low dust- 
emission coefficient (DEC, the total emission of the dust to the 
total emission of the star in the wavelength range 7.7-14.0/im) 
for S Ori of 0.24, i.e. similar to RR Sco (DEC=0.21) or R Leo 
(DEC=0.23), indicating an optically thin dust shell. 

3. VLTI/MIDI measurements 

We obtained mid-infrared (8-1 3 //m) interferometry of S Ori 
with the instrument MIDI (Leinert et al. 2003) at the ESO VLT 
Interferometer (VLTI; Glindemann et al. 2003J in service mode 
between 31 December 2004 and 30 December 2005. Light was 
combined from two 8 m Unit Telescopes (UTs). We chose to 
use the prism with spectral resolution R - A/AA = 30 to 
disperse the interferograms. Two various techniques of beam 
combinations were used, namely HIGH_SENS (HS, 2004-12-31 
to 2005-03-03) and SCI.PHOT (SP, 2005-1 1-12 to 2005-12-30) 
combinations. Using the HS combination, all arriving light is 
first combined interferometrically, and the photometric spec- 
trum for each incoming beam is recorded sequentially. Using 
the SP combination, beam splitters are used to record the in- 
terferograms and the photometric spectra simultaneously. An 
estimate of the photometric count rates at the time of record- 
ing the interferometric data is needed to compute the raw vis- 
ibility values, and the SP technique is thus expected to result 
in a higher precision of the visibility measurement for bright 
targets. Our observations followed the standard procedures for 
MIDI service mode observations as described in detail in the 
MIDI user manuaQ. For principles of observations with MIDI 
see also Leinert et al. (120041) . 

The details of our observations are listed in Table [T] Listed 
are the epoch, date, time, and Julian Date (JD) of the observa- 
tion; the visual phase Oyis; the baseline configuration used with 
its ground length B; the dispersive element and beam combiner 
BC (HS for High_Sens and SP for Sci_Phot); and the pro- 
jected baseline length Bp and its position angle PA on sky (deg. 
east of north). The last two columns describe the ambient con- 
ditions DIMM seeing, and coherence time tq, both at 500 nm. 
Dates of observation were combined into 4 epochs A, B, C, 
D. Each epoch has a width of less than 17 days, i.e. less than 
5% of S Ori's variability period. The dates of observation and 
the epochs compared to the visual lightcurve are indicated in 
Fig.[T] We estimate the uncertainty in assigning a visual phase 
value to each of our observing epochs to ~ 0.1, taking into ac- 
count both the time widths of our epochs and the uncertainty in 
the choice of S Ori's current variability period and date of last 
maximum (see Sect. |2]l. Epochs A and B are located around 
stellar minimum with a separation of about 2 months (phase 
difference ~0.1), while epochs C and D occurred about a year 
later shortly after the following stellar maximum with a sep- 
aration of 6 weeks. Additional MIDI/UT data on S Ori from 
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Table 1. VLTI/MIDI observations of S Ori. 



Ep. 


Date 


Time 


JD 




Conf. 


B 


Disp. 


BC 


B„ 


PA 


Seeing 


To 






UTC 


-2.4E6 








Elem. 




[m] 


[deg] 


["] 


[msec] 


A 


2004-12-31 


02:43 


53370.6 


0.42 


UT3-UT4 


62 m 


Prism 


HS 


62.27 


-71.9 


0.76 


3.4 


B 


2005-02-21 


02:37 


53422.6 


0.54 


UT2-UT4 


89 m 


Prism 


HS 


74.31 


83.0 


0.74 


3.6 


B 


2005-03-01 


02:39 


53430.6 


0.56 


UT3-UT4 


62 m 


Prism 


HS 


41.65 


-55.6 


0.62 


3.6 


B 


2005-03-03 


00:59 


53432.5 


0.56 


UT2-UT4 


89 m 


Prism 


HS 


84.35 


82.7 


0.98 


3.7 


C 


2005-11-12 


08:57 


53686.9 


1.16 


UT1-UT4 


130 m 


Prism 


SP 


123.02 


62.8 


0.54 


4.4 


C 


2005-11-15 


07:21 


53689.8 


1.16 


UT2-UT3 


47 m 


Prism 


SP 


46.26 


43.9 


0.61 


2.9 


D 


2005-12-30 


03:00 


53734.6 


1.27 


UT1-UT4 


130 m 


Prism 


SP 


123.89 


59.4 


1.12 


2.8 



Table 2. Properties of the MIDI calibration stars used. 



HD 


Name 


Sp.T 


Old 




-f 12 


<r 








[mas] 


[K] 


[Jy] 


[°] 


50778 


eCMa 


K4III 


3.90±0.02 


4049 


24.6 


27 


31421 


0^ Ori 


K2III 


2.65±0.01 


4500 


13.6 


12 


49161 


17 Men 


K4III 


2.44±0.01 


4248 


7.2 


20 


61935 


a Men 


G9 III 


2.24+0.01 


4813 


10.2 


36 


94510 


u Car 


Kim 


2.23+0.01 


4951 


6.9 


90 


107446 


eCru 


K3.5III 


4.43±0.02 


4139 


32.4 


101 



2004-12-29, 2005-01-02, 2005-02-28, 2006-03-12, 2006-03- 
12, as well as MIDI/AT data taken between October 2005 and 
March 2006 were not used because of various technical prob- 
lems with these data sets. The interferometric field of view of 
MIDI used with the UTs is ~ 250 mas, which corresponds to 
about 25 times the expected photospheric angular diameter of 
SOri. 

We used CMa (HD 50778) as the interferometric calibra- 
tor, a K4 giant at a distance on sky of 27°, which was observed 
within about 30 minutes of each S Ori observation. A number 
of observations of other MIDI calibration stars using the same 
instrument settings were obtained during some of our observa- 
tion nights for other programs, and we make use of these data 
as well. Table |2] lists the calibration stars and their properties. 
These parameters are taken from the instrument consortium's 
catalog^. Listed are the HD number, the name, the spectral type, 
the limb-darkened diameter ©ld, the effective temperature T^s, 
the flux at 12jum, and the distance on sky from S Ori ^. The 
UD diameter for K giants in the mid-infrared based on spheri- 
cal model atmospheres is smaller by ~ 5% than ©ld (estimated 
with procedures as used in Wittkowski et al. '2006'). HD 50778 
is our main calibrator The other calibrators were observed dur- 
ing some of our observation nights for other programs. Among 
the interferometric calibration stars listed in Table|2l absolutely 
calibrated spectra are available for HD 31421 and HD 61935 
(Cohen et al. [T999b . 

3.1. MIDI data reduction 

Mean coherence factors and photometric count rates were 
obtained for each data set of S Ori and the calibration 
stars using the MIA+EWS software package, version 1 .4 (Jaffe, 

^ http://www.ster.kuleuven.ac.be/~tijl/MIDI_calibration/mcc.txt 



Koehler, Cotton, Chesneau, et al0), complemented by the tool 
MyMidiGui, version 1.41 (Hummell). Data for 32 spectral 
channels within the nominal MIDI wavelength range from 
8.0jum to 13.0yum were used. The detector masks were calcu- 
lated by the procedure of MyMidiGui, and were used for both 
the MIA and the EWS analysis. The results obtained by MIA and 
EWS were found to be virtually identical for all datasets. We 
chose to use in the following the results based on the EWS anal- 
ysis. It is based on a coherent integration of the fringe signal af- 
ter numerical compensation of optical path differences in each 
scan. 

Calibrated visibility spectra for S Ori were derived using 
the transfer function obtained for our calibrator HD 50778. The 
scatter of all obtained transfer functions during the night was 
used to estimate its uncertainty. The latter depends on the sam- 
ple of calibrations stars available per night, and we assumed 
a minimum uncertainty of the transfer function of 0.1 based 
on nights when many values were obtained. The total error of 
our calibrated S Ori visibility spectrum includes the errors of 
the raw coherence factors computed by EWS, the uncertainty of 
the calibrator diameters, and the scatter of the transfer function 
during the night. The total error is dominated by systematic 
uncertainties for the height of the interferometric transfer func- 
tion, so that the relative accuracy of the visibility values among 
the spectral channels is better than indicated by our total error 
bars. 

We absolutely calibrated the photometric spectrum of S Ori 
for each data set using the interferometric calibrator HD 50778 
also as the photometric calibration star The absolutely cali- 
brated photometric spectrum of HD 50778 was not available 
directly. We used instead the average of the available spectra 
of HD 20644 and HD 87837 (Cohen et al. |T999l ). scaled with 
the IRAS 12;um flux (Beichmann et al. 119881 ) to the level of 
HD 50778. These stars have the same spectral type (K4III) and 
very similar effective temperatures as HD 50778 (4058 K and 
4094 K versus 4049 K). In addition, we calibrated the photo- 
metric spectrum of HD 50778 using the spectro-photometric 
calibrator HD 31421, which was observed close in time during 
the night 2005-03-03, and verified that our synthetic HD 50778 
spectrum is valid. We obtained integrated 8.0-13.0/zm S Ori 
flux densities for our MIDI epochs A to D of 187 ± 12 Jy, 
152 ± 17 Jy, 198 + 18 Jy, and 203 + 13 Jy. For comparison, 



^ http://www.strw.leidenuniv.nl/~koehler/MIDI 
http://www.sc.eso.org/~chummel/midi/midi.html 
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the integrated IRAS LRS flux density for the same bandpass 
(8.0-13.0;um) is 169 Jy. 

The MIDI calibrated visibility, as well as calibrated flux 
values, are detailed in Table ??. 

3.2. MIDI data reduction results 

As an alternative and probably more intuitive quantity com- 
pared to the visibility values, we calculated the equivalent UD 
diameters. The equivalent UD diameter was computed for each 
projected baseline length and spectral channel separately and 
defined as the diameter of that uniform disc that results in the 
same visibility value at each separate data point. Note that the 
true intensity profile is not expected to be a uniform disc. As 
a result, the equivalent UD diameters for different projected 
baseline lengths are not expected to coincide (not even for a 
circularly symmetric intensity profile). However, these values 
are fully equivalent to the calibrated visibility values and can 
provide a first rough estimate of the characteristic size of the 
target. 

Figures |2]-(5] show, for MIDI epochs A-D, respectively, the 
obtained MIDI flux spectra, the MIDI visibility values as a 
function of wavelength, the MIDI visibility values as a func- 
tion of spatial frequency for the example of three spectral chan- 
nels and the equivalent UD diameter values. Also shown are the 
synthetic visibility and flux values of the atmosphere and dust 
models described below in Sect. 13.31 in each panel, as well as 
the intensity profiles of these models for the example of three 
bandpasses in the bottom panels. 

The general shape of our S Ori visibility functions, or 
equivalent UD diameter values, in Figs. [JHH is qualitatively 
similar to the MIDI data of the Mira variable RR Sco described 
in Ohnaka et al. (120051 ). Objects such as these are character- 
ized by a partially resolved stellar disc, including atmospheric 
molecular layers that are optically thick at mid-infrared wave- 
lengths, surrounded by a spatially resolved, optically thin dust 
shell. The visibility increases between about 8-9 ;um, corre- 
sponding to a quasi-constant equivalent UD diameter that is 
roughly twice the continuum photospheric size. Here, the ob- 
served intensity is dominated by radiation from optically thick 
molecular layers at distances up to roughly two continuum pho- 
tospheric radii. Beyond about 9-I0//m, the visibility function 
flattens, corresponding to an increasing equivalent UD diame- 
ter. In this wavelength region, spatially resolved, optically thin 
radiation from the dust shell starts to represent a considerable 
part of the total intensity. In addition, extinction of the stellar 
light by the dust shell becomes important. At the longest wave- 
lengths of MIDI, spatially resolved radiation from the dust shell 
dominates the measured intensity. A detailed model of S Ori's 
atmosphere and dust shell is described in Sects. [331 and l3.4l 

3.3. MIDI data modeling 

Mid-infrared photometric and interferometric data of Mira 
stars are in general sensitive to details of the structure of the 
stellar atmosphere in which molecular layers of various geo- 
metrical and optical characteristics (e.g., Tej et al. I2003al) lie 



above the continuum-forming photosphere, as well as to the 
surrounding dust shell (cf. e.g., Ohnaka et al. 120051 |2006bt 
Ireland & Scholz |2006l Weiner et al. |2006] l. We used the P and 
M model atmosphere series (Ireland et al. 2004a, 2 004bl ) to 
describe the stellar atmosphere including the continuum pho- 
tosphere and overlying molecular layers, and added an ad-hoc 
radiative transfer model of the dust shell. The details of our 
modeling are described in the following. 

3.3.1 . Atmosphere model for S Ori 

Few dynamic atmosphere models for oxygen-rich Mira stars 
are available that include the effects of molecular layers (see, 
e.g., Scholz l200n 12003) . The P and M model series (Ireland 
et al. I2004al I2004bl l are complete self-excited dynamic model 
atmosphere of Mira stars, and they have successfully explained 
a lot of observational data (e.g.; Hofmann et al. ,2000. 20011 
l2002t Woodruff^ et al. l2004t Fedele et al. 12005]) . These model 
series have been constructed to match the prototype oxygen- 
rich Mira stars o Cet and R Leo. The P and M series differ 
with respect to the mass of the so-called "parent star", which is 
the hypothetical non-pulsating equivalent of the pulsating Mira 
variable. The geometric pulsation of the Mira occurs around the 
parent star's Rosseland radius R - ^(tross = 1) (Ireland et al. 
I2004all2004bl l. The parent star has solar metallicity, luminosity 
L/Lo=3470, period 332 days, mass M/Mq^I.O (P series) and 
1.2 (M), radius Rp/RQ^24l (P) and 260 (M) (Hofmann et al. 
119981 ). The moderately larger radius of the M-series parent star 
leads to a slightly lower effective temperature for the parent 
star and systematically lower phase-dependent effective tem- 
peratures for the pulsating Mira compared to the P series. The 
M models tend to exhibit less pronounced cycle-to-cycle varia- 
tions than the P models and to have more compact atmospheres 
(Ireland et al. l2004bl l. Effects on interferometric diameter mea- 
surements resulting from the differences between the P and M 
model series are subtle. When comparing P and M models of 
similar phases to near-infrared interferometric data of R Leo, 
the resulting Rosseland angular diameter values were found to 
agree within < L5% (Fedele et al. l2005] l. 

Compared to o Cet and R Leo, for which the P and M 
model series are designed, S Ori is a slightly cooler Mira 
variable (M6.5-M9.5 versus M5-M9/M6-M9.5, Samus et al. 
120041) with a longer period (420 d versus 332 d/3 10 d, Samus et 
al. 120041), a larger main sequence precursor mass (1.3 Mq ver- 
sus 1.0 Mo/I .2 Mo, Wyatt & Cahn [T983T l, and a larger radius 
(mean continuum photospheric radii roughly R ~450Ro ver- 
sus -350 Ro; based on BW05, Woodruff^ et al. l2004l Fedele et 
al. 120051 ). However, when scaled to variability phases between 
and 1 and to the corresponding angular size on sky, the gen- 
eral model results are not expected to be dramatically different 
for S Ori compared to Mira stars, such as o Cet and R Leo. 
Hereby, the M series is expected to be better-suited to describ- 
ing S Ori than the P series because of its higher parent-star 
mass. We concluded that the use of the M series is currently 
the best available option for describing S Ori's atmosphere in- 
cluding the continuum photosphere and overlying atmospheric 
molecular layers. 
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Fig. 2. Results of our VLTI/MIDI 8- 13 pm interferometry of S Ori at 
epoch A, stellar phase 0.42. Shown are our MIDI data compared to 
the best-fitting model with parameters listed in Tabled In each panel, 
the x-symbols with error bars denote the measured values, the solid 
lines the synthetic values of the overall model consisting of the stellar 
model atmosphere and the dust shell. The dotted lines denote the con- 
tribution from the attenuated atmosphere model alone, and the dashed 
lines the contribution of only the dust shell (attenuated input radiation 
subtracted). 

a) The calibrated MIDI flux spectrum. 

b) The calibrated MIDI visibility values as a function of wavelength. 

c) The calibrated MIDI visibility values as a function of spatial fre- 
quency for the example of three spectral channels at wavelengths 
8.5 yum, 10.5 yum, and 12.5 fim. 

d) The equivalent UD diameter values. Also indicated are the pho- 
tospheric angular diameter Qphot and the inner dust shell boundary 
©i„ as obtained by our model fit, as well as the mean 42.8 GHz and 
43.1 GHz SiO ring angular diameters (Sect.Q. 

e) The model intensity profiles at spectral bands 8-9 yum, 1 0-1 1 fim, 
and 12-13 fim. The bottom abscissa shows the radius as a function of 
the photospheric radius ^phot, and the top abscissa the corresponding 
angular scale obtained from our model fit. The arrows indicate the lo- 
cation of the inner boundary radius of the dust shell as well as of 
the mean 43. 1 GHz and 42.8 GHz SiO maser spot radii i?43.i and Rn.i- 
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Fig. 3. Results of our VLTI/MIDI 8— 13/ini interferometry of Fig. 4. Results of our VLTI/MIDI 8— 13//ni interferometry of 
S Ori at epoch B, stellar phase 0.55. As Fig. [2] S Ori at epoch C, stellar phase L16. As Fig. [J] 
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Fig. 5. Results of our VLTI/MIDI 8-13yum interferometry of 
S Ori at epoch D, stellar phase 1.27. As Fig.|2] 



3.3.2. Dust-shell model for S Ori 

Except for the dust-enshrouded very-long-period Mira IW Hya 
(Jeong et al. 120031) . models that self-consistently combine de- 
tailed dynamic atmosphere models and the formation of dust 
shells are currently not readily available for oxygen-rich AGB 
stars. However, the first attempts in this direction have just been 
made (Ireland & Scholz'2006', Woitke'2006^. 

Dust shells around AGB stars have often been modeled 
using ad-hoc radiative transfer calculations (see; for instance; 
Danchi et al. '1994'; Lorenz-Martins & Pompeia'2000; Ohnaka 
et al. ,2005. ,2006a), and we followed this approach as the only 
currently available option for describing S Ori's dust shell. We 
employed the Monte Carlo radiative transfer code mcsimjnpi 
(Ohnaka et al. 2006a) to calculate the dust shell model of S Ori. 
For the central radiation source of the radiative transfer model, 
we used the spectral energy distribution (SED) of the respective 
M model atmosphere for the wavelength range from 0.32 /urn 
to 23yum and extended it by a blackbody approximation for 
longer wavelengths. Dust temperatures were calculated in ra- 
diative equilibrium. 

Our MIDI spectra of S Ori do not exhibit any prominent sil- 
icate feature, indicating that silicate grains are not a major con- 
stituent of the dust shell. Lorenz-Martins & Pompeia ( I2000I I 
have investigated the envelopes of 3 1 oxygen-rich AGB stars 
using the IRAS LRS spectra. They classify the objects into 
three groups, whose IRAS LRS spectra can be modeled using 
either silicate grains, AI2O3 grains, or a mix thereof. S Ori is 
part of their study and belongs to the group of objects that can 
be modeled well with AI2O3 grains alone. They show that a 
blend of alumina AI2O3 grains (Koike et aI. [T995T l for A < Sfxm 
and porous amorphous AI2O3 grains (Begemann et al. 1 19971 1 
for A > 8fim gives the best agreement with the IRAS spec- 
tra, so we followed this approach (as in Ohnaka et al. 12005 1 for 
the case of RR Sco). Likewise, we considered small amounts 
of warm silicates (Ossenkopf et al. ,1992) following the same 
works. The grain size was set to 0.1 /zm for both AI2O3 and 
silicates. We used the same inner boundary radius and density 
distribution for both dust species. The amount of dust of each 
species was characterized by the optical depths at a wavelength 
of 0.55 ;um, denoted Tv(Al203) and Tv(silicate). The density 
distribution p(r) was described by a single power law with in- 
dex p ipif) oc r^P). We used the inner boundary radii in units 
of the stellar photospheric radius Rphot- All radiative transfer 
models of the dust shell were calculated with an outer radius, 
so that the shell thickness is RoatlRin - 1000. 

3.3.3. The main model parameters 

The stellar atmospheric structure including the effects from at- 
mospheric molecular layers is fully described by the choice of 
the specific model of the M series. The M series consists of 20 
models for different phase and cycle combinations as detailed 
in Ireland et al. (I2004bl l. Here, we used a grid of the six M mod- 
els of cycle 2 of the series, which are M20 (model visual phase 
0.05), M21n (0.10), M22 (0.25), M23n (0.30), M24n (0.40), 
and M25 (0.50). These models cover the range of our obser- 
vations at phases 0.16 (epoch C), 0.27 (epoch D), 0.42 (epoch 
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A), and 0.55 (epoch D). The absolute model phases with re- 
spect to observations are uncertain by at least 0.1, with relative 
uncertainties of 0.01 to 0.02 between different models of the 
series (Ireland et al. l2004b] l. Likewise, the absolute assignment 
of visual phases to our epochs of observation is uncertain by 
about 0. 1 (see Sect.O. For each of these 6 atmosphere models, 
whose SEDs are used as central source, we calculated a grid of 
radiative transfer models of the dust shell. These grids included 
all combinations of optical depths Ti/(Al2O3)=1.0, 1.5, 2.0, 2.5, 
3.0; ri/(silicate)=0%, 10% of Ti/(Al203); inner boundary radii 
i?in/^Phot=l-8, 2.0, 2.2, 2.4, 2.6; and gradients /:'=2.0, 2.5, 3.0, 
3.5,4.0. 

Any overall model was thus defined with this set of five in- 
dependent parameters: (1) choice of the dust-free M model, (2) 
TviAhOi), (3) Tv(siHcate), (4) Ri^, (5) p. We used the well- 
defined continuum photospheric radius at A =1.04jum Rim 
of the respective dust-free M model to characterize the pho- 
tospheric radius /Jphot; i-S- we set /Jphot = ^i.04- Any overall 
model is rigid; i.e. the distance of any model layer from ^phot is 
fixed with the parameters given above. The corresponding an- 
gular diameter ©phot was treated as the only fit parameter when 
comparing any such overall model to our MIDI data. 

3.3.4. Synthetic visibility and flux values 

Stellar contribution Monochromatic CLVs of the M models 
were calculated at 3000 equally-spaced frequencies for the 
wavelength range 8//m< A < 13/zm as described in Ohnaka 
et al. (I2006bl l. covering radii r between and R^^^^ - 5 x Rp. 
These CLVs were calculated with H2O and SiO (28SiO, 29SiO, 
and 30SiO) lines included. The line lists used for the calcula- 
tions were the same as those used for modeling RR Sco by 
Ohnaka et al. (W05); and as in the case of RR Sco, we adopted 
siHcon isotope ratios of 28Si/29Si = 15 and 28Si/30Si = 20. 
The H2O and SiO line opacities were computed assuming local 
thermodynamic al equilibrium (LTE) and a Voigt profile with a 
micro-turbulent velocity of 3 km s"' . 

The visibility values based on the stellar contribution were 
calculated for each spectral channel c of MIDI as 



rj^f' ir(r) e-r^^h dr dv. 



(1) 



Here, B is the projected baseline length used for the observa- 
tion. Note that this computation is valid for the case of the EWS 
data reduction, which coherently integrates the fringe signal 
to directly obtain the visibility modulus. For the MIA analysis, 
which is a power spectrum analysis integrating the squared vis- 
ibility modulus, the integration over frequency would need to 
be performed over the squared monochromatic visibility values 
(cf. Wittkowski et al.'2006'). The observed stellar flux contribu- 
tion f^^'' was computed for each of MIDI'S spectral channels 



as 



^star _ o ^ 

Tc - 2 7r 



Qout["]/2 
206265 



^v,„i„(c) Jo 



Cir) e-^^'W dr dv. (2) 



The spectral channels of MIDI were assumed to have a rect- 
angular sensitivity function between frequencies Vmin and Vmax- 



The term e describes the extinction by the dust shell. The 
optical depth t(v) was calculated based on tabulated absorption 
and scattering coefficients Qabs and Qsca of the two dust species 
/ as 



t(v) = > r(Vo) : : , 

where we used vo corresponding to Aq - 0.55/im. 



(3) 



Dust contribution Using the dust temperatures and monochro- 
matic mean intensities computed by a Monte Carlo run with 
mcsiirunpi for a given set of dust shell parameters, we cal- 
culated monochromatic CLVs of the dust shell (excluding the 
contribution of the central star) with the ray-tracing method 
as described in Ohnaka et al. (I2006al I2007I I. The radii r cov- 
ered the range between and 257?phot, and we denote the 
outermost model radius, where we set the intensity to 0, by 
^ouf = 25 X /?phot- The outer radius of 25 /?phot was chosen be- 
cause the MIDI interferometric field of view of ~ 250 mas cor- 
responds to 25 times the expected photospheric angular diam- 
eter of S Ori of ~ 10 mas. We verified that changing the outer 
model radius from 25/?phot to 507?phot changed the diameter 
results by less than 0.3%. This indicates that the contribution 
from dust to the overall CLV outside of 25 photospheric radii 
can be neglected and that the exact choice of the outermost 
model radius is not critical for our calculation. Visibility (V^"^') 
and flux (Z^"**') values for each spectral channel of MIDI were 
computed in the same way as for the dust-free model atmo- 
sphere, but without the extinction term e"^*''\ 

Comparison of visibility and flux values to our data Finally, 
the flux and visibility values of the overall model were com- 
puted as 



fr-fr+fc 



•dust 



and 



total 



Y. 



■dust 



(4) 



(5) 



ytotal y*total 

This simple addition of the CLVs of the dust-free model at- 
mosphere of the M series and of the pure dust-shell (i.e. after 
subtraction of the contribution by the attenuated input radia- 
tion) is fully valid as long as the stellar atmosphere and dust 
shell are spatially separated. In fact, there may be some over- 
lap of the extended low-intensity atmospheric molecular layers 
of the dust-free atmosphere model and the dust shell in our 
overall models. However, we verified that the inner boundary 
radius of the dust shell is clearly larger than that of the layer at 
which the atmosphere becomes optically thick at mid-infrared 
wavelengths. Ireland & Scholz (I2006I I find for similar condi- 
tions that the gas temperatures for a dust-free M model and an 
M model that includes re-heating of gas by dust, differ by less 
than 50 K for radii up to about 3 parent-star radii. 

For any given overall model separately, we optimized the 
angular diameter ©phot using a standard least-square fit that 
minimizes the value between synthetic and observed visi- 
bility and flux values. Hereby, ©phot was treated as the only 
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Table 3. Fit results of our MIDI measurements of S Ori. 



Epoch (<I)vis) 


A (0.42) 


B (0.55) 


C(1.16) 


D(1.27) 


Best model 




M22 


M24n 


M23n 


M21n 


•tModel 




0.25 


0.40 


0.30 


0.10 


0Model - 0Vis 




-0.17 


-0.15 


-1-0.14 


-0.17 


Best RJRphot 






1.8 


2.0 


2.2 




2.4 


Best TvCAlzOa) 




2.5 


2.5 


1.5 




1.5 


Best Tv/(Silicate) 




















Best p 






3.5 


3.5 


3.0 




2.5 


©Phot [mas] 




9.00 


7.92 


9.66 


9.48 






0.16 


0.10 


0.43 


0.35 


n 






62 


186 


124 




62 


cTformai [mas] 




0.10 


0.05 


0.02 


0.07 


o"modd [mas] 




0.30 


0.09 


0.06 


0.34 


Table 4. Reduced xl 


values obtained for our 


grid of over 


models. 


















Epoch 






A 


B 


C 




D 




M model 


M20 




0.27 


0.16 


0.52 




0.38 






M21n 




0.26 


0.14 


0.53 




0.35 






M22 




0.16 


0.17 


0.81 




0.51 






M23n 




0.17 


0.12 


0.43 




0.61 






M24n 




0.20 


0.10 


0.66 


1.04 






M25n 




0.19 


0.17 


0.53 




0.58 




^in/^Phot 


1.8 




0.16 


0.10 


0.56 




0.54 






2.0 




0.19 


0.10 


0.52 




0.46 






2.2 




0.21 


0.11 


0.43 




0.39 






2.4 




0.24 


0.12 


0.50 




0.35 






2.6 




0.29 


0.18 


0.57 




0.38 




Tv'CAlzOa) 


1.0 




/ 


/ 


0.77 




0.48 






1.5 




0.18 


0.11 


0.43 




0.35 






2.0 




0.17 


0.11 


0.56 




0.54 






2.5 




0.16 




0.10 


0.80 


0.71 






3.0 




0.18 


0.11 


/ 


/ 


Tv(Silicate)/ 


0% 




0.16 




0.10 




0.43 




0.35 




Tl.(Al203) 


10% 




0.19 


0.12 


0.52 




0.38 




p 


2.0 




0.31 


0.20 


0.87 




0.38 






2.5 




0.20 


0.14 


0.54 




0.35 






3.0 




0.17 


0.10 


0.43 




0.39 






3.5 




0.16 




0.10 


0.50 


0.51 






4.0 




0.20 


0.12 


/ 


/ 



free parameter. The resulting reduced;i'^ values among our full 
grid of different overall models were compared in order to find 
the best-fitting set of model parameters. This analysis was sep- 
arately performed for each of our epochs A-D. 

3.4. MIDI model results 

Table |3] lists for each of our IVIIDI epochs the parameter set of 
the best-fitting overall model of our grid. Listed are for each of 
our 4 epochs the best-fitting atmosphere model of the M series; 
its model phase; the difference between observed and model 
phase; the best-fitting parameters of the dust shell (tv(A1203), 
Tv(Silicate), Tv(Silicate), p); the best-fitting angular diameter 
©Phot; the reduced;^y value; the number of data points used (n); 



and the formal and model errors. The model errors are derived 
as the standard deviations of ©phot values based on different 
overall models of our grid that are consistent at the 1 cr level 
based on application of the F-test. The model errors of epochs 
A and D are larger than those of epochs B and C, which can be 
explained by the use of only one baseline configuration at these 
epochs. Table|4]lists for each epoch and for each grid parameter 
the best xl value that could be obtained by any combination 
of all other parameters. The best xl value for each epoch is 
marked by boxes. 

The synthetic visibility and flux values corresponding to 
the best-fitting overall models are shown in Figs. |2]|5]in com- 
parison to the measured values. Panels (d), which show the 
equivalent UD diameter values, also indicate the photospheric 
angular diameter 0phot and the inner dust shell boundary ©in 
derived from the model fit. Panels (e) show the intensity pro- 
files of the overall models that fit our MIDI data best for three 
bandpasses at 8-9 yum, 10-1 1 yum, and 12-13yum. These panels 
have a bottom and a top abscissa indicating (bottom) the radius 
as a function of the model photospheric radius ^phot and (top) 
the corresponding angular scale using the best-fitting angular 
diameter ©phot- Panels (e) also indicate the inner dust shell ra- 
dius Rin, as well as the mean SiO maser ring radii R42.S. and 
/?43.i, which are derived in Sect.|4] A two-dimensional pseudo- 
color image of the mid-infrared model intensity showing the 
molecular layers and the dust shell is displayed in Fig. |2l The 
photospheric disk is represented by a light blue color Overlaid 
are the images of the SiO maser radiation described in Sect. HI 

The photospheric diameter ©phot The obtained photospheric 
diameter values ©phot at phases 0.42, 0.55, 1.16, and 1.27 are 
9.0 mas + 0.3 mas, 7.9 mas + 0.1 mas, 9.6 mas + 0.1 mas, and 
9.5 mas + 0.4 mas, respectively. Here, the error bars include the 
formal, as well as the model, errors from Table|3] Our measure- 
ments indicate significant changes of the photospheric diameter 
©Phot as a function of stellar phase. A detailed comparison of 
©Phot to the inner dust shell radius and the SiO maser ring radii 
follows in Sect. |5] 

In BW05, we estimated S Ori's photospheric angular di- 
ameter at phase 0.73 to ~ 9.2 mas, based on our /T-band UD 
diameter measurements and using correction factors from UD 
to photospheric diameters from Ireland et al. (I2004al I2004bl l. 
This value is consistent with the photospheric diameters de- 
rived here, as the estimate of 9.2 mas at phase 0.73 from BW05 
lies as expected in between the values of 7.9 mas at phase 0.54 
and 9.6 mas at phase 1.16 derived in the present paper 

The /T-band UD diameter of 10.54 mas + 0.68 mas by van 
Belle et al. ( 119961 ) was obtained at phase 0.56. At phases close 
to this value, the M series results in a ratio between the A'-band 
UD diameter and the photospheric diameter of ~ 1 .35 (Fig. 7 of 
Ireland et al. |2004b] l. With this correction, the photospheric di- 
ameter corresponding to this measurement at phase 0.56 would 
be 7.8 mas + 0.5 mas, and would thus be fully consistent with 
our value of 7.9 mas + 0. 1 mas at the same phase (Epoch B at 
phase 0.54). 

MiUan-Gabet et al. ( 120051 ) obtained a /T-band UD diame- 
ter of 9.6 mas + 0.2 mas at phase 0.1. At this phase, the ra- 
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tio between /T-band UD diameter and photospheric diameter is 
predicted to be close to unity (Ireland et al. |2004bT l. This mea- 
surement thus also agrees well with our photospheric diameter 
of 9.6 mas + 0.4 mas at phase 1.16 (epoch C). 

The good agreement of the photospheric angular diameters 
derived from our MIDI measurements in the present paper with 
previous near-infrared diameter measurements increases confi- 
dence in the validity of our modeling using the M model series 
complemented by a radiative transfer model of the dust shell. 

The choice of the M model The intensity profile of the ex- 
tended atmosphere, and thereby the choice of the best-fitting 
model out of the M series, is mainly constrained by the shape of 
the visibility functions, in particular for longer baseline lengths, 
in the wavelength range ~ 8-9 fj.m, where the contribution by 
the dust shell is relatively low. 

The differences in visual stellar phase of the best-fitting M 
model compared to the visual phase at the epochs of obser- 
vation are -0.17, -0.15, -1-0.14, and -0.17, respectively. These 
differences may be due to the uncertainties of the absolute as- 
signments of visual phases to the model series as well as to 
the dates of observation, which are both on the order of 0.1. 
The relative phase values are expected to agree better, which is 
true for epochs A, B, and D. Differences between these phases 
might also occur owing to the different radius and pulsation pe- 
riod of the model series' parent star compared to S Ori. These 
parameters might have an effect on the extension of the at- 
mosphere, so that the best-fitting extension of the atmosphere 
might belong to a model at a different visual phase. An effect 
such as this might explain the clearly different phase difference 
between observation and model at epoch C compared to the 
three other epochs. 

The good agreement of the shape of the synthetic visibility 
function with our observations at the wavelength range from 
8 fim to 9 jum and the simultaneous agreement with the mea- 
sured flux level again confirms that the M model series is usable 
for S Ori. 

As an indication of the extension of the molecular layers 
of the best-fitting M models, the 50%/10% intensity radii of 
theA^-band (8-13yum) model CLVs are 1.5/1.6, 1.8/1.9, 1.6/1,7, 
and 1.4/1.9 photospheric radii Rptot for epochs A-D, i.e. for 
models M22, M24n, M23n, and M21n, respectively. 

The dust shell chemistry The ratio of the number of AI2O3 
grains to that of silicate grains is mostly constrained by the 
shape of the flux and visibility spectra at the silicate feature 
close to IQfim. We obtain best fits to our S Ori data with a dust 
shell only consisting of AI2O3 grains. The addition of small 
amounts of silicate grains does not improve the fits for any of 
our epochs. While our grid in Table |4] only includes silicate 
grains with the same inner boundary radius as for the AI2O3 
grains, we also compared to our data the best-fitting models 
from Table[3]with silicate grains added at larger inner boundary 
radii of 4 and 6 photospheric radii. These models lead to clearly 
increased values as well. In summary, our data indicate that 
silicon is not bound in dust, neither at very close distances of 
~ 2 Rphot nor at larger distances of 4-6 Rphot- A dust shell only 



consisting of AI2O3 grains is consistent with the modeling of 
the IRAS LRS spectra of S Ori by Lorenz-Martins & Pompeia 
(12000b . 

The radial density profile of the dust shell The total opti- 
cal depth of the dust shell is mainly constrained by the height 
and shape of the MIDI total flux spectrum. The radial density 
profile described by the inner boundary radius and the den- 
sity gradient is mostly constrained by the shape of the visi- 
bility function versus wavelength. Here, the values at shorter 
wavelengths, probing warmer dust, constrain the inner bound- 
ary radius best. The density gradient is mostly constrained by 
the visibility values at longer wavelengths and shorter baseline 
lengths, with which the dust shell is less resolved. 

We obtain best-fitting inner boundary radii of 1.8, 2.0, 2.2, 
and 2.4 photospheric radii /?phot at our phases 0.42, 0.55, 1.16, 
and 1.27, respectively. These radii are well-constrained and do 
not significantly change for variations in other model param- 
eters. We estimate the uncertainty of this parameter to about 
0.2 X 7?phot. The formation of dust consisting of AI2O3 grains at 
these short distances from the stellar surface is consistent with 
the empirical results by Lorenz-Martins & Pompeia ( 20001 1. as 
well as with the recent theoretical calculations by Ireland & 
Scholz (I2OO6I 1 and Woitke (l2006l l. 

Best-fitting density gradients p at the same phases are 3.5, 
3.5, 3.0, and 2.5, respectively, and best-fitting optical depths 
Ty are 2.5, 2.5, 1.5, and 1.5. Taking the inner boundary radii, 
the density gradients, and the optical depth values together, our 
fits indicate more compact dust shells with larger optical depth 
near the stellar minimum and more extended dust shells with 
lower optical depth after stellar maximum. 

Wind models We used the radiation transport code DUSTY 
(Ivezic & Elitzur 119971 Ivezic et al. 11999 ) with its option 
of modeling radiatively driven winds to compute the density 
structure of S Ori's envelope that would result from a wind 
model. We computed wind models for each of our epochs A 
to D using a blackbody with effective temperature of the best- 
fitting M model as the central radiation source and with the Ri^ 
and Ty values from Table [3] We used the same optical proper- 
ties of AI2O3 as above, and we set Rom/Rin to 10^. Figure |6] 
shows the resulting normalized density profiles of the wind 
models together with the p oc r^^ and p oc r^^ curves. The wind 
models result in density structures p oc r^^ for r ^ 10/?phot. 
Closer to the photosphere, the gradients of the wind density 
profiles become steeper than p oc r"^. As our MIDI visibility 
measurements are mostly sensitive to the inner dust region, the 
relatively high density gradients given in Table [3] could be ex- 
plained by such a wind model. The differences among the em- 
pirically found density gradients for the different epochs could 
be explained by different mean distances from which most mid- 
infrared photons escape. The wind models shown correspond 
to mass-loss rates between ~ 2.4 lO"*" Mg/yr (Epoch C) and 
~3.5 lO^^'Mg/yr (Epoch B). The indication of more compact 
dust shells with larger optical depth near stellar minimum and 
more extended dust shells with lower optical depth after stellar 
maximum thus suggests a wind with increased mass loss and 
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Table 5. VLBA observations of S Ori. 
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Fig. 6. Density profiles of radiatively driven wind models cal- 
culated with the DUSTY code (Ivezic & Elitzur 119971 Ivezic 
et al. 1 19991 1 for input parameters from Table |3] together with 
p oc r'^ and p oc r^^ curves. 



dust formation occurring relatively close to the stellar surface 
near minimum visual phase and with an expanded dust shell 
after maximum phase. 

Symmetry of the molecular and dust shells Our MIDI data at 
epochs A and D are each obtained at one single position angle 
on sky. Epoch B spans a range in position angles of ~ 40 ° and 
epoch C of ~ 20 °. Within these probed ranges, our S Ori data 
can be described well by spherically symmetric distributions of 
molecular layers and dust shells without the need to introduce 
deviations from circular symmetry to our modeling. However, 
we cannot rule out asymmetries at wider ranges of position an- 
gles, in which case the resulting parameters in Table[3]are only 
vahd along the position angles on the sky given in Table [T] 
Note that the distributions of the SiO maser spots toward S Ori 
shown in BW05, as well as below in Sect. 21 do not indicate a 
globally asymmetric gas distribution. 

4. VLBA measurements 

Contemporaneous with the VLTI/MIDI observations, we ob- 
served the y = 1,7= 1-0 (43.1 GHz) and v = 2, 7 = 1 - 
(42.8 GHz) SiO maser emission associated with S Ori (a = 
05''29"'00^9,5 =-04°41'32".7, J2000) using the 10 stations of 
the VLBA, which is operated by the National Radio Astronomy 
Observatory (NRAO). Table|5]lists the details of our SiO maser 
observations recorded over three epochs. Listed are the epoch, 
date, time. Mean Julian Date (MJD) of the start of observation, 
experiment duration, number of antennas, and the visual phase 
Ovis- The first two epochs occurred near stellar minimum and 



Ep. 


Date 


Time 


MJD 


Dur. 


No. 


*vis 






UTC 


-2.4E6 


(hrs) 


Ant. 




A 


2005-01-17 


02:36 


53387.6 


5.0 


9 


0.46 


B 


2005-02-28 


23:47 


53430.5 


2.6 


10 


0.56 


C 


2005-11-05 


07:23 


53679.8 


5.0 


10 


1.14 



were spaced 42 days apart or approximately 10% of the stellar 
period. The third epoch occurred roughly 9 months later just 
after maximum on the next stellar cycle. Data were recorded in 
a similar manner for all three epochs at reference frequencies 
of 43.122080GHz and 42.820587GHz for the v = 1 and v = 2 
SiO transitions, respectively. Right- and left-hand ckcular po- 
larization was used with 8-MHz (56.1 km s ') bands centered 
on the local standard-of-rest (LSR) velocity of 18.0km s^^ 
Interspersed with observations of S Ori were 10-minute scans 
of extragalactic continuum sources (0423-013 and 0359-1-509) 
for the purpose of delay and bandpass calibration. For epoch 
A only 9 of the 10 stations of the VLBA were available for 
our observations. Unfortunately, the one antenna that was un- 
available was the station at Mauna Kea, HI. This station when 
combined with the St. Croix, VI antenna provides the longest 
baseUne (greatest resolution) in the array. For this reason the 
resolution of the epoch A images is reduced from that of the 
other two epochs. There was also a problem with the second 
epoch (epoch B) of observations in that roughly half of the 
data were accidentally lost after the correlation. Despite this 
loss, we were still able to use the remaining data to produce 
quality images of the masers. NRAO kindly scheduled a third 
epoch of observations (epoch C) to make up for the loss. 

4.1. VLBA data reduction and analysis 

The data were correlated at the VLBA correlator operated by 
NRAO in Socorro, New Mexico. Auto and cross-correlation 
spectra consisting of 256 channels with channel spacings 
of 31.25kHz (~0.2kms"') were produced by the correla- 
tor Calibration was performed using the Astronomical Image 
Processing System (AIPS) maintained by NRAO. The data 
were calibrated in a manner similar to what was performed in 
BW05. For each epoch, residual delays due to the instrumenta- 
tion were corrected by performing a fringe fit on the continuum 
calibrator scans. Residual group delays for each antenna were 
determined and applied to the spectral line data. Variations in 
the residual delays ranged from 2-4 ns resulting in phase errors 
of no more than 1.5-3° across the 8-MHz band. 

The bandpass response was determined from scans on the 
continuum calibrators and used to correct the target-source 
data. The time-dependent gains of all antennas relative to a ref- 
erence antenna were determined by fitting a total-power spec- 
trum (from the reference antenna with the target source at a 
high elevation) to the total power spectrum of each antenna. 
The absolute flux density scale was established by scaling these 
gains by the system temperature and gain of the reference an- 
tenna. Errors in the gain and pointing of the reference antenna 
and atmospheric opacity variations contribute to the error in 
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the absolute amplitude calibration, which is accurate to about 
15-20%. 

Residual fringe rates were obtained by fringe-fitting a 
strong reference feature in the spectrum of each maser tran- 
sition. For epoch A we used the channel at velocity Vlsr = 
12.3 km s ' for the v = 1 transition and the channel at velocity 
Vlsr = 16.9 km s ' for the v - 2 transition. For epoch B we 
used channels at Vlsr = 12.8kms"' and Vlsr - 12.1 kms"' 
for the V = 1 and v - 2 transitions, respectively. Finally, for 
epoch C we used the same velocity channel for both transi- 
tions Vlsr = 10.3 kms '. The resulting fringe-rate solutions 
were applied to all channels in each spectrum. An iterative self- 
calibration and imaging procedure was then performed to map 
this reference channel for each transition. The resulting resid- 
ual phase and amplitude corrections from the reference chan- 
nels at 42.8 GHz and 43.1GHz were applied to all channels 
in the respective bands. This procedure was repeated for each 
transition at each epoch. 

In order to accurately compare the distributions of the two 
maser transitions, it is desirable to determine a common spa- 
tial reference point. However, after the fringe-fitting step to de- 
termine residual fringe-rates, all absolute position information 
is lost for the VLBA data. To accomplish this registration, we 
used the same technique employed in BW05 with the residual 
fringe-rates, phases, and amplitudes from one transition (e.g., 

V - 2, 42.8 GHz) used to calibrate the other transition (e.g., 

V - 1, 43.1 GHz) in order to determine an offset between the 
two phase centers. 

The final image cubes of the SiO maser emission consisting 
of 1024 X 1024 pixels (~51 x 51 mas) were generated for each 
transition at each of the three epochs. Images were produced 
for all spectral channels containing SiO maser emission in each 
transition at each epoch. Synthesized beam sizes used for epoch 
A were 0.59 x 0.40 mas and 0.60 x 0.37 mas for the v = 1 and 

V — 2 transitions, respectively. For epoch B these sizes were 
0.52x0.21 mas and 0.5 1x0.20 mas. For epoch C the beam sizes 
used in the imaging were 0.46 x 0.16 mas and 0.44 x 0.16 mas. 
Note that the resolution was greatly reduced in epoch A due 
to the loss of the Mauna Kea, HI VLBA antenna mentioned 
previously. The resolution of epoch B is slightly worse than 
that of epoch C due to the loss of half the data, thus reducing 
the uv coverage of the array. Using the procedure described 
above, the image cubes from the two transitions were spatially 
aligned for each epoch. Off-source rms noise estimates in the 
images ranged from 5 mJy to 9 mJy for epoch A, 3 mJy to 
7 mJy for epoch B, and 7 mJy to 18 mJy for epoch C. Figure|7] 
shows the total intensity images of the v = 1, 43.1 -GHz and 

V = 2, 42.8-GHz SiO maser emission toward S Ori overlaid 
onto a pseudo-color representation of the mid-infrared intensity 
profile obtained in Sect. [3] 

The six resulting image cubes (2 transitions at 3 epochs) 
were then analyzed to extract relevant maser parameters. Two- 
dimensional Gaussian functions were fit to maser emission 
above a cut-off flux density of 100 mJy in each spectral (ve- 
locity) channel using the AIPS task SAD. This fitting yielded 
velocities, flux densities, and positions in right ascension and 
declination for all emission components identified in the im- 
ages. Errors in right ascension and declination of identified 



components were computed using the fitted source size divided 
by twice the signal-to-noise ratio (s/n) in the image and ranged 
from 1 yuas for features with high s/n, to 28 yuas for features with 
lower s/n. 

The remaining analysis of the maser component identifica- 
tions was performed outside of the AIPS package. For com- 
parison with the total flux densities derived from the MIDI 
data, we summed the flux densities and their errors for all fit- 
ted components in all channels for each transition/epoch. These 
integrated flux densities are discussed further in Sect. |5] Since 
the -0.2 kms ' channel spacing is sufficient to resolve the 
masers spectrally, features typically appear in multiple adjacent 
spectral channels. Positions in right ascension and declination 
and center velocities for the masers were determined using a 
flux-density-squared weighted average for features identified 
in two or more adjacent channels with a spatial coincidence 
of 0.5 mas. The flux assigned to the maser averages was the 
maximum single-channel flux density. The results from the im- 
age analysis and component averaging are described in the next 
section. 

4.2. VLBA results 

4.2.1 . The spatial structure of the SiO masers 

The maser features characterized using the procedures de- 
scribed above are listed in Tables ??-?? and are shown in 
Fig. |8] where there are six panels ordered by epoch left to right 
and by transition top to bottom. Within each panel there are 
two sub-panels. Each upper sub-panel shows the spectrum of 
the maser emission with flux density plotted as a function of 
local standard of rest (LSR) velocity. Point sizes are propor- 
tional to the logarithm of the flux density and are color-coded 
according to velocity bin in increments of 1.7kms The bot- 
tom sub-panels of Fig. [8] show the spatial distribution of the 
SiO masers as observed in the images. Point sizes are again 
proportional to the logarithm of the flux density with the same 
velocity color-coding as in the top sub-panels. 

Examining the upper sub-panels and the component tables, 
we find that the masers have shifted slightly toward the blue 
end of the spectrum from our previous BW05 observations. 
There the maser velocities ranged from 11.1 to 24.7 kms"', 
while for these latest three epochs the masers span a velocity 
range from 8.6 to 19.0 km s"' . The velocity range for our latest 
three epochs is consistent with earlier spectral measurements 
by Jewell et al. (I1991I I. who measured a velocity range of 7.1 
to 19.6 km s~' for the v = 1 masers. Spectral ranges measured 
by Cho et al. (1 19961 1 were much narrower at 1 1. 0-16. 0km s ' 
for the V = 1 and 12.0-14.7kms ' for the v - 2 masers, re- 
spectively. For the three epochs discussed here, flux densities 
of the peak component in the spectrum ranged from 3.0 to 
5.4Jybeam ' for the v = 1 masers and 1.3 to 6.4Jybeam~' 
for the V - 2 masers. These are close to our BW05 results for 
which we measured peaks of 4.3 Jy beam ' and 5.3 Jy beam"' 
for the V = 1 and v = 2 masers, respectively. 

When compared to the spectrally-averaged total-intensity 
images (Fig.|7]i, the maser component maps in lower sub-panels 
of Fig.[8]provide an accurate representation of the spatial struc- 
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ture of the maser emission summed over all velocity channels 
in the image cube. The component maps all show that the SiO 
masers in both transitions form partial to full rings of emission 
with the typical clumpy distribution. The top three panels show 
the V - 1,7 = 1-0 to form a partial ring for epoch A with 
a gap to the north, nearly a full ring for epoch B , and again 
a partial ring dominated by emission in the northern and west- 
em regions for epoch C. For the v - 2,J - 1-0 SiO masers 
shown in the bottom three panels, we observe nearly a full ring 
of emission for epoch A, a partial ring with a gap to the south 
for epoch B, and a partial ring dominated by northern and west- 
ern features for epoch C. In all three epochs the v = 1 masers 
appear to be distributed over a larger area on the sky than the 
V = 2 masers. In addition, the v = 1 rings appear thicker than 
the corresponding v - 2 rings at each epoch. It should be noted 
that, while epochs A and B were spaced only 42 days apart, 
there are significant differences in the number and distribution 
of the masers in both transitions. This is not unexpected since 
previous multi-epoch VLBl studies (e.g. Diamond & Kemball 
I2003I I have shown the fine structure of the SiO masers to vary 
significantly over periods as short as two weeks. This variation, 
combined with the relative complexity of the maser structures, 
precluded us from tracking individual maser components in or- 
der to determine maser proper motions. 

To better characterize the apparent size and thickness of 
the SiO maser distributions, we computed the average angu- 
lar distance of the masers from the center of their distribution 
{?). To accomplish this, we first determined the center of the 
distribution by performing a least-square fit of a circle to the 
combined v = 1 and v - 2 maser positions. We were able to 
combine the two sets of masers only because we had previously 
aligned the images of the two transitions following the proce- 
dure described in Sect. 14.11 This fit produced a common center 
from which we computed the mean maser angular distance rsio 
and the standard deviation for each transition independently for 
each epoch. 

The mean angular distances from the center for the ob- 
served SiO masers at v = 1, 43.1 GHz and v = 2, 42.8GHz 
are listed in Table |6] These distances are also indicated by a 
dashed circle in each of the six panels in Fig. [8] The standard 
deviations in the mean distances are also listed in Tableland 
represented in Fig. [8] as the error on the mean. Two times the 
standard deviation also provides an indication of the thickness 
of the shell. The listed parameters in Tableware for each epoch 
of observation/visual phase and each of the two transitions, the 
apparent mean angular distances from the center of the masers 
spots r,,=i,2, the ratio of the mean distances of the two transi- 
tions r,,=i/r,,=2, as well as the mean LOS velocities Vy=\;i and 
their ratio Vv=i/Vy=2- 

From Table |6] we see that ? ranges from 8.6 to 10.2 mas. 
At every epoch we find that the v = 1 masers have a greater 
mean distance from center than the v = 2 masers. This find- 
ing is consistent with BW05 in which we found r,,=i -9.4 mas 
and r,,=2 - 8.8 mas. In BW05, however, we were unable to 
make a definitive statement about the relative sizes of the two 
maser distributions because we had only the one epoch of ob- 
servations, and at this epoch, the v = 2 masers were primarily 
confined to a small region on the northwest side of the shell. In 



Table 6. Results from our VLBA observations of the v = 1, 
43.1 GHz and v = 2, 42.8 GHz SiO maser transition toward 
S Ori. 



Epoch 




A 




B 


C 


Phase Ovis 




0.46 




0.56 


1.14 


[mas] 


9.9 


± 1.4 


9.6 


± 1.1 


10.2 ± 1.2 


ry^2 [mas] 


8.6 


± 1.0 


9.1 


± 1.1 


8.9 ± 0.9 






1.15 




1.05 


1.15 


V,^i [kms-'] 


12.6 


±2.4 


13.3 


±2.5 


13.0 ±2.8 


^,,=2 [kms-'] 


12.5 


±2.4 


12.9 


±3.2 


12.5 ± 2.9 






1.01 




1.03 


1.04 



our more recent epochs, the rings of emission are much more 
complete making the comparison between transitions easier. 
Although the differences in ? for epochs A, B, and C are less 
than the standard deviation (thickness of the ring), they are 
still greater than the resolution of the images given earlier by 
the synthesized beam sizes of 0.6 x 0.4 mas for epoch A and 
0.5 X 0.2 mas for epochs B and C. The fact that the v = 1 ring 
is consistently larger than the v - 2 ring suggests that for S Ori 
the two sets of masers are not co-spatial. 

The relative separation of the v = 1 and v - 2 masers for 

5 Ori is also consistent with recent results obtained for the SiO 
masers around other stars including: IRC H- 1001 1 (Desmurs et 
al. l2000] l: TX Cam (Desmurs et al. 20001 Yi et al. l2U05T l: o Cet, 
U Ori, and R Aqr (Cotton et al. 2006). In addition, numeri- 
cal simulations of SiO masers (Humphreys et al. 119961 Gray 

6 Humphreys 2000 ) show the v = 1,7 - 1-0 masers to 
occur farther from the star than the v = 2,7=l-0 masers 
throughout the entire pulsation. The simulations of Gray & 
Humphreys (I2000I I also show that the mean thickness of the 
V = 1 shell is roughly twice that of the v - 2 shell throughout 
the stellar cycle. In our measurements, we find the v = 2 shell 
to be thicker in the BW05 observations, the v = 1 shell to be 
thicker in epochs A and C, and the two transitions to be roughly 
equivalent in epoch B. In the two cases for which the v = 1 
shell was thicker, the ratio of the two thicknesses is 1.3-1.4, 
much less than the 2 determined by Gray & Humphreys. The 
angular sizes of the maser shells are reported in Table|2l At the 
assumed distance for S Ori of 480 pc + 120 pc (van Belle et al. 
120021 ). 1 mas is roughly equivalent to 0.48 AU. Thus the linear- 
scale sizes for maser distributions range from 8.3 to 9.8 AU and 
the maser shell thickness (twice the standard deviation) ranges 
from roughly 0.9 to 1.3 AU. 

4.2.2. The kinematics of tine SiO masers 

In BW05, we reported that the SiO masers showed no coher- 
ent velocity structure indicative of global expansion/infall or 
rotation. Examining Fig.[8l we again find no evidence of such 
coherence in the line-of-sight (LOS) velocities of the masers as 
a function of spatial location. We did, however, notice that there 
appears to be a velocity gradient at all epochs, with masers 
toward the blue- and red-shifted ends of the spectrum lying 
closer to the center of the distribution than masers at intermedi- 
ate velocities. This phenomenon was also found in BW05 and 
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Fig. 9. Histogram of the LOS velocities of our SiO maser spots 
combining VLBA epochs A-C. Panel (a) shows the histogram 
for the V - 2 42.8 GHz masers and panel (b) for the v = 1 
43.1 GHz masers. Also indicated are the mean values y,,=i and 
y,,=2- Vv^i is slightly redder than Vv=2- 

has been previously observed in other SiO maser studies (e.g. 
Boboltz & Marvel 2000! Hollis et al. 200TT). 

The LSR stellar velocity of S Ori has not been measured 
particularly well. In the General Catalog of Stellar Radial 
Velocities (Wilson il953) . the value listed for S Ori is 22 km s ' . 
Young (I1995I I observed the sub-millimeter CO (3-2) and CO 
(4-3) lines toward S Ori and determined values of Vlsr = 
14.5 + 0.2 and 14.1 + 0.5 km s"' from fits to the spectra, re- 
spectively. Winters et al. ( I2003I I derived Vlsr - 14.0 km s"' 
based on fits of the CO (1-0) and CO (2-1) lines. Since the SiO 
masers are not distributed evenly about any of the above veloc- 
ities, we computed an average LOS velocity from the masers 
themselves for each transition at epoch. The values of Vlos for 
each of the epochs are reported in Table |6] with errors of 15% 
estimated from the uncertainty of the absolute amplitude cali- 
bration. Although the differences are well within the error bars, 
it appears that the v - 2 masers are shifted slightly toward the 
blue end of the spectrum relative to the v = 1 masers at all three 
epochs. In addition to computing values for individual transi- 
tions and epochs, we also combined the three epochs and com- 
puted a mean LOS velocity for each transition. We obtained 
mean values of Vlsr = 12.6 + 2.8 kms"' and 13.0 + 2.6 km s"' 
for the 42.8 GHz and 43.1 GHZ masers, respectively. These 
values are shown in Fig. |9] Again we see that Vi^i is slightly 
redder than Vy=2- 

To verify the notion that the higher-velocity masers are 
found closer to the star, we plotted component distance from 
the center of the distribution (radius) versus its velocity relative 
to the average LOS velocity for that transition (Vlos - V'los)- 
The results are shown in Figs.[TO]and[TT]for the v - l,J- 1-0 



43.1GHz and the v = 2,7 = 1 -0 42.8 GHz transitions, re- 
spectively. In the figures we plotted each epoch with a differ- 
ent symbol. The masers in both transitions appear to have the 
same distribution with a central peak near Vlos - Vlos = and 
decreasing maser radius as a function of increasing values of 
I Vlos - Vlos I- 

A simple model that is often used in the case of OH (e.g. 
Reid et al. rT977] l and H2O (e.g. Yates & Cohen [T994l l maser 
kinematics is that of a uniformly expanding thin shell. In this 
model the projected distance from center (a) of a maser on the 
shell is related to its LOS velocity V^ - V* by the expression: 



'"sh 



v-v* 

Vsh 



= 1, 



(6) 



where r^h and Vsh are the radius and the expansion velocity of 
the shell, respectively. The LSR velocity of the star (V*) is 
assumed here to be equal to Vlos- Numerical simulations of 
SiO masers often employ the Sobolev or large velocity gradi- 
ent (LVG) approximation (Lockett & Elitzur .1992. Bujarrabal 
[19941 Doel et al. [19951 Humphreys et al [2002] l. The LVG ap- 
proximation assumes a large velocity gradient across the SiO 
maser region, thus producing large Doppler shifts that serve to 
disconnect a region of velocity coherence from the surround- 
ing medium (Bujarrabal 1994). The parameter typically used 
to characterize the velocity field in the maser region is the log- 
arithmic velocity gradient e given by 



d\nV 



e = 



r dV 



(7) 



d\nr 

A value of e = is equivalent to constant velocity expansion, 
while e > 1 corresponds to a velocity field with large radial ac- 
celerations. The SiO maser numerical simulations (e.g. Lockett 
& Elitzur il992, Bujarrabal , 1994. Doel et al. il995; ) typically use 
a value of e = 1 . Chapman & Cohen d 19861 1 determined values 
of e for the various maser species surrounding the supergiant 
VX Sgr and found e 1 in the SiO region, e » 0.5 in the 
region of the H2O and mainline (1665 MHz, 1667 MHz) OH 
masers, and e a; 0.2 in the 1612 MHz OH maser region. 

Plotted along with the data in Figs. [10] and [TTj are four 
curves computed using the above equations. They represent 
shells at the inner and outer boundaries of the SiO masers. 
For the masers in each transition, we used an inner shell ra- 
dius equivalent to the smallest average distance minus its stan- 
dard deviation from Table|6] Similarly, for the outer shell radius 
we used the largest average distance plus its standard deviation 
from Table [6] For the scenario of a logarithmic velocity gradi- 
ent, the velocity at the inner v = 2 boundary was set at 7 km s"' 
and the velocities of the other shells computed assuming e = 1 . 
The velocity of the outermost shell of the v = 1 transition is 
Vsh ~ 10.5kms-i. 

For the scenarion of a constant velocity expansion (e = 0), 
the best by-eye fit results in Vsh - 10km s"'. In this case, the 
higher expansion velocity at the inner boundary of the v = 2 
transition is required in order to provide a reasonable fit to the 
outer boundaries of both transitions. For the e - I case, the 
inner v - 2 velocity can be lower with the radial acceleration 
allowing for a better fit to the outer boundaries. For compari- 
son, the escape velocities computed for the innermost v = 2 and 
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Fig. 10. Distances of the 43.1 GHz maser components from the 
common center of the distribution (radius) versus their velocity 
relative to the average LOS velocity for the 43.1 GHz transi- 
tion (Vlos - Vlos)- Maser components of the different epochs 
are denoted by different symbols and colors. The plot indicates 
that higher-velocity masers are found closer to the star. This 
trend can be explained by means of expanding maser shells. 
The curves represent shells at the inner and outer boundaries 
of the SiO maser shell using (solid lines) a logarithmic veloc- 
ity gradient and (dashed lines) constant velocity expansion. See 
text for details. 



outermost v = 1 shells, assuming a mass for S Ori of 1.0 Mq, 
are 22.1 and IS.Okms"', respectively. Therefore, even at the 
velocities bracketed by the curves, the maser gas is still gravi- 
tationally bound. 

From the SiO maser data on S Ori, we cannot definitively 
rule out constant velocity expansion; however, the SiO maser 
region is very dynamic, influenced by the pulsation of the star 
and the passage of shocks. Theoretical studies all seem to in- 
clude some form of velocity gradient across the SiO maser re- 
gion. In addition, constant velocity expansion would be more 
conducive to masers along the line of sight to the star because 
of the increased velocity coherence. SiO masers are rarely ob- 
served in front of the stellar disc, and for S Ori in particular, 
such masers have not be seen. We therefore consider some form 
of acceleration in the region the more likely scenario. 

The canonical interpretation of the rings formed by the SiO 
masers is that the masers are tangentially amplified and primar- 
ily confined to a region close to the plane of the sky. Under this 
interpretation, the ring size and thickness are direct indicators 
of the equivalent 3-D shell size and thickness. The radially ex- 
panding spherical shell described above will naturally result in 
a projection effect where masers out of the plane of the sky have 
smaller angular distances from center due to the projection. In 
the case of S Ori, this projection is rather small. The equation 



- (km s-') 



Fig. 11. As Fig. [To] but for the 42.8 GHz transition. The axes' 
scales are identical for the two figures. Compared to Fig. [TOl 
it is evident that the 42.8 GHz maser spots lie at systematically 
closer distances to the center of the common maser distribution 
than the 43.1 GHz maser spots. 

for the angle 9 measured between the LOS and a point on the 
shell locating the maser can be written in terms of the velocities 

as 

1/2 

(8) 



sin 6 - 



1 - 



V_ 



In the most extreme case, all of the thickness observed in the 
ring is due to projection of the masers on the outermost thin 
shell with a velocity of Vsh ~ 10.5kms"' (e - 1). Figures|9]- 
[TTIshow that the majority of the masers are within ~4kms"^ 
of LOS center velocity. For these extreme 4 km s"' masers, 
sinO ~ 0.92 and 6 » 67°. This results in a diff'erence of ~8% 
between the projected distance of the maser and its true radial 
distance. For masers less than 4 km s"' the difference is smaller. 
Assuming the masers are distributed uniformly, the differences 
between the projected ring sizes and thicknesses and their 3- 
D equivalents is ~4%. Thus, the given errors of the projected 
ring sizes determined from the maser distributions include the 
uncertainty of the projection effect. 

5. VLTI/MIDI and VLBA results as a function of 
stellar phase 

Table [T] provides an overview on the photometric and spatial 
parameters of S Ori as a function of observational epoch/stellar 
phase, as derived from both our VLTl/MlDl and VLBA obser- 
vations in Sects. [3] and |4] The listed photometric parameters 
include the V magnitude based on the lightcurve in Fig.lT] the 
A^-band flux based on our photometric MIDI data (Sect. 13. Il l, 
and the total flux values of the two maser transitions (Sect. l4Tt . 
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Table 7. Overview of the photometric and spatial parameters 
of S Ori as a function of stellar phase, derived from both our 
VLTl and VLBA observations. 





A 


B 


C 




D 


Phase Ovis 


0.44 


0.55 


1.15 




1.27 


V [mag] 


12.8 ±0.3 


14.0 ± 0.3 


9.4 ± 0.3 


10.5 


±0.4 


In [Jy] 


187 ± 12 


152 ± 17 


198 ± 18 


203 ± 13 


743.1 [jyj 






iqc + Qf) 

lyo m j\j 






/42.8 [Jy] 


148 ± 22 


46 ± 7 


383 ± 57 






©Phot [mas] 


9.0 ± 0.3 


7.9 ±0.1 


9.7 ±0.1 


9.5 


± 0.4 


011, [mas] 


16.2 ± 1.9 


15.9 ± 1.6 


21.2 ± 1.9 


22.8 


±2.1 


043.1 [mas] 


19.6 ±2.8 


19.0 ± 2.6 


20.4 ± 2.4 






042.8 [mas] 


19.0 ±2.0 


17.8 ± 3.4 


17.8 ± 1.8 






0iii/0Phot 


1.8 ±0.2 


2.0 ± 0.2 


2.2 ± 0.2 


2.4 


±0.2 


043.1 /0Pliot 


2.2 ± 0.3 


2.4 ± 0.3 


2.1 ±0.3 






042.8 /0Phot 


2.1 ±0.2 


2.3 ± 0.4 


1.9 ± 0.2 






043.1 /0i„ 


1.2 ±0.2 


1.2 ±0.2 


1.0 ± 0.1 






042.8/0111 


1.2 ±0.2 


1.1 ±0.2 


0.8 ± 0.1 






^Phot [Ro] 


469 ± 127 


411 ± 111 


498 ± 134 


491 : 


t 133 



The listed spatial information includes the photospheric angu- 
lar diameter Ophot (Sect. I3.4l i. the inner dust shell diameter ©!„ 
(Sect. l3.4T i. and the maser ring angular diameters ©43.1 and ©42.8 
(Sect. 14.2. 11 1. To compare the continuum photospheric diame- 
ter, the inner dust shell diameter, and the maser shell diame- 
ters, the ratios ©in/Ophot, ©43.i/42.8/©Phot, and ©43. 1/42.8 /©in ai'e 
listed as well. As an absolute scale, the stellar photospheric ra- 
dius /?phot at each phase is derived from ©phot and the adopted 
distance to S Ori (Sect. |2]). Here, the variability phases are the 
mean values of the respective phases of the IVIIDI and VLBA 
observations. 

The continuum photospheric radii Rphot, the inner dust shell 
radii Rin (corresponding to ©in) and the maser shell radii 7^43 1 
and /?42.8 (corresponding to ©43.1/42.8) are indicated in panels 
(d) of Figs. 12115] which show the synthetic mid-infrared CLVs. 
This allows a comparison of these radii to each other and to the 
profile of the extended molecular atmosphere. Figure |7] pro- 
vides as well a graphical comparison of the photospheric disc, 
the molecular layers, the dust shell, and the maser spots for 
each epoch. 

Variability of visual, mid-infrared, and maser flux densities 
The visual, mid-infrared, and maser fluxes all show clear vari- 
ability as a function of observational epoch. The stellar phases 
given throughout this paper are visual phases as introduced in 
Sect. |2] The A^-band flux shows clear variability with an 
amplitude of about 50 Jy (~ 30%). /n exhibits a minimum at 
our epoch B at visual phase 0.55, i.e. shortly beyond the visual 
minimum, rises toward stellar maximum, and continues to in- 
crease between epochs C and D at visual phases 1.15 and 1.27, 
probably indicating a phase lag and different shape of the A^- 
band variability curve compared to the visual lightcurve. The 
total maser fluxes are significantly lower near stellar minimum 
(epochs A and B) compared to stellar maximum (epoch C), 



with the exception of the strong /43 1 value at epoch B. The lat- 
ter can most likely be explained by the larger number of maser 
spots forming a nearly complete ring-like structure for this tran- 
sition and epoch, while the other transitions and epochs show 
only partial maser rings. 

Photospheric and dust shell radii as a function of visual 
phase The photospheric angular diameter ©phot derived from 
our MIDI observations and modeling shows a clearly phase- 
dependent size with an amplitude of about 1 .8 mas (~ 20%) that 
is well-correlated in phase with the visual lightcurve. It shows a 
minimum at epoch B, which is nearest to the visual minimum, 
and a maximum at epoch C, which is nearest to the visual max- 
imum. Also the inner dust-shell radius shows a variation (am- 
plitude ~ 7 mas ~ 30%), which is more closely correlated in 
phase and amplitude to than to V. The inner dust shell ra- 
dius is smallest at our epoch B near visual minimum and largest 
at epoch D beyond visual maximum (phase 1.27). Compared 
to ©Phot, the inner dust shell boundary is located significantly 
closer to the stellar surface at our epochs A and B near visual 
minimum (1.8-2.0x©phot) than at the post-maximum visual 
phases C and D (2.2-2.4 x©phot)- Simultaneously, as discussed 
in Sect. 13.41 the density gradient p is significantly steeper and 
the optical depth significantly higher near visual minimum than 
at the post-maximum visual phases. This can be understood by 
a wind model with increased mass loss near minimum phase 
occurring close to the stellar surface and an expanded dust shell 
after maximum (cf. Sect. 13. 4I) . 

Maser ring radii as a function of visual phase The mean SiO 
maser ring radii obtained at our three VLBA epochs show a 
small variability with amplitude of < 7%, which is only a frac- 
tion of the mean width of the maser shells of ~ 13% and clearly 
below the variability amplitude of the photospheric and dust 
shell radii at the same epochs. However, we cannot rule out 
that the maser ring radii show a greater variability in between 
our epochs of observation. Observations by Cotton et al. (I2006l l 
consistently find variations in maser ring diameters of 3-14% 
for different Mira stars, but no clear correlation with visual vari- 
ability phase. Theoretical estimates by Humphreys et al. ( I2002I I 
predict a variability of the 43. 1 GHz maser shell radius with an 
amphtude of ~ 20%. 

Maser ring radii compared to the photospheric radius The 

true location of the star relative to the maser spots is unknown. 
Based on the assumption that the center of the maser spot distri- 
bution coincides with the center of the star, a comparison of the 
maser ring radii with the continuum photospheric radii gives 
an estimate of the distance of the maser spots from the stel- 
lar photosphere. The ratios between the maser ring radii and 
the photospheric radii in Table [7] range between 2.1+0.3 and 
2.4+0.3. These values are consistent with the values reported in 
BW05 (1.9 and 2.0), with those by Cotton et al. (l2004b . which 
range for different Mira stars -not always compared at the same 
phase- between 1.8 and 2.9, and the value by Fedele et al. 
(120051) for R Leo of 1.8. Our observations indicate a variabiUty 
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in the ratio of the mean maser shell radius to the photospheric 
radius with an amplitude of ~ 10% for the 43.1 GHz transition, 
and ~2Q% for the 42.8 GHz transition, with uncertainties on 
the same order This observed variability is dominated by the 
variability of the stellar continuum photospheric diameter. 

Maser ring radii compared to tine molecular atmosphere 
Panels (d) of Figs. |2}{5]as well as Fig. |7]show that the mean 
maser ring radii at all our epochs mark the region just beyond 
the steepest decrease in the mid-infrared model intensity, i.e. 
shortly outward of the densest region of the molecular atmo- 
sphere. 

Maser ring radii compared to the inner dust shell bound- 
ary Table |7] as well as panels (d) of Figs. llHll indicate that 
the inner boundary of the AI2O3 dust grains near visual min- 
imum (epochs A and B) appears to be located between the 
steepest decrease in the mid-infrared model intensity and the 
mean location of the SiO maser spots. This means that the 
low-density atmospheric molecular layers, the SiO maser spots, 
and the AI2O3 dust grains are co-located near visual minimum. 
At our post-maximum epochs C and D, the inner boundary of 
the AI2O3 dust shell appears significantly expanded, while the 
maser shells remain at about the same location. As a result, 
the inner boundary of the AI2O3 dust grains appears to have 
expanded outward of the mean SiO maser rings at the post- 
maximum phase. 

6. Summary, conclusions, and discussion 

We have obtained the first time series of observations of a 
Mira variable star that include concurrent infrared and radio in- 
terferometry. For the present paper, we obtained mid-infrared 
interferometry of S Ori with VLTl/MIDI at four epochs in 
December 2004, February/March 2005, November 2005, and 
December 2005. We observed v = 1 , 7 = 1 - (43. 1 GHz) and 
V - 2,J - 1-0 (42.8 GHz) SiO maser emission toward S Ori 
with the VLBA in January, February, and November 2005, 
concurrent to within 5% of the variability period to the first 
three VLTI/MIDI epochs. The first two common VLTl/VLBA 
epochs are located near visual minimum (phases 0.44 and 0.55) 
and the third after next visual maximum (phase 1.15). The 
fourth MIDI epoch occurred at a later post-maximum phase 
of 1.27. 

The mid-infrared MIDI data are sensitive to the structure of 
the atmosphere consisting of the continuum photosphere and 
overlying molecular layers, as well as to the properties of the 
dust shell. The MIDI visibilities and flux densities are modeled 
using the recent M model atmosphere series by Ireland et al. 
j2004bl l. which is a dust-free self-excited dynamic model at- 
mosphere series, to which we add an ad-hoc radiative transfer 
calculation of the dust shell. The constrained model parame- 
ters include the continuum photospheric angular diameter, as 
well as characteristics of the dust shell. The MIDI visibility 
and total flux data at all epochs can be described well by the 
chosen approach of combining the dust-free M model atmo- 
sphere series with an ad-hoc radiative transfer model of the 



dust shell. Diff'erent model parameters can be constrained by 
diff'erent subsets of our MIDI data. 

The resulting continuum photospheric angular diameters 
at visual variability phases 0.42, 0.55, 1.16, and 1.27 are 
9.0 + 0.3 mas, 7.9 + 0.1 mas, 9.7 + 0.1 mas, and 9.5 + 0.4 mas, 
respectively. The photospheric radius thus shows a significantly 
phase-dependent size with amplitude of ~ 20% that is well cor- 
related in phase with the visual lightcurve. 

The dust shell can be modeled well with AI2O3 grains 
alone, and our data show no indication of silicates, consis- 
tent with the modeling of the IRAS LRS spectra by Lorenz- 
Martins & Pompeia (I2000I I. The inner dust shell angular di- 
ameters are 16.2 ± 1.9mas, 15.9 + 1.6 mas, 21.2 ± 1.9mas, and 
22.8 ±2.1 mas, respectively, for the visual phases given above, 
corresponding to 1.8, 2.0, 2.2, and 2.4 times the continuum 
photospheric diameters. Thus, dust starts to form very close 
to the stellar surface, corresponding to the region of the low- 
intensity wings of the atmospheric molecular layers of our 
dust-free atmosphere model. In particular, the inner dust shell 
radius is smallest at the near-minimum phases, located directly 
outside of the steep decrease in the dust-free, mid-infrared 
model intensity. The inner dust-shell boundary expands to 
larger distances from the stellar surface at our post-maximum 
phases. The overall AI2O3 optical depth, as well as the den- 
sity gradient, is largest at our near-minimum phases and signif- 
icantly lower at the post-maximum stellar phases. These dust 
shell characteristics can be explained by a radiatively driven 
wind with increased mass-loss near visual stellar minimum oc- 
curring close to the stellar surface, and an expanded dust shell 
at the post-maximum phases. 

The 43.1GHz and 42.8 GHz maser spots show the typi- 
cal structure of partial to full rings with a clumpy distribu- 
tion. There is no sign of a globally asymmetric gas distribution. 
The 43.1GHz mean maser ring angular diameters at visual 
phases 0.46, 0.56, and 1.14 are 19.6 ±2.8 mas, 19.0 ±2.6 mas, 
and 20.4 ±2.4 mas, respectively. The values for the 42.8 GHz 
maser are systematically smaller by 3-19% and amount to 
19.0 ±2.0 mas, 17.8 ±3.4 mas, and 17.8 ± 1.8 mas. The indica- 
tion of systematically smaller mean maser ring diameters of 
the 42.8 GHz transition compared to the 43.1GHz transition 
is consistent with earlier results in the literature on other Mira 
stars. The velocity structure of both maser transitions indicate 
an expansion of the maser shell. This expansion is most likely 
an accelerated expansion, but we cannot completely rule out a 
uniform expansion. The expansion velocity for both transitions 
ranges between about 7kms"' and about lOkms ' between 
the inner and outer shell radii based on a scenario of an ac- 
celerated expansion or amounts to about 10km s"' based on a 
scenario with uniform expansion. 

Compared to the continuum photospheric radius /?phot, 
the 43.1 GHz/42.8 GHz SiO maser ring radii are located at 
2.2 ±0.3/2.1 ±0.2 (phase 0.44), 2.4 ±0.3/2.3 ±0.4 (0.55), and 
2.1 ±0.3/1.9±0.2 (1.15) times /?phot- These values are consis- 
tent with theoretical estimates and earlier observations avail- 
able in the literature. Compared to the molecular atmosphere, 
the maser shells mark the region of the molecular atmosphere 
just beyond the steepest decrease in the mid-infrared model 
intensity. Compared to the dust shell, the maser spots are co- 
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located with the inner dust shell near stellar minimum. At our 
post-maximum phases, the maser spots remain at about the 
same location, while the dust shell has expanded outward. 

Altogether, our observations indicate a pulsation of the 
stellar photosphere approximately in phase with the visual 
lightcurve. Our measurements suggest an increased mass-loss 
rate near stellar minimum and formation of AI2O3 dust grains 
occurring close to the stellar surface at about 1.8-2.0 photo- 
spheric radii directly outward of the dense molecular atmo- 
sphere, co-located with the extended low-intensity wings of the 
molecular atmosphere and with the SiO maser spots. At the 
post-maximum phases, the dust shell has expanded to about 
2.2-2.4 photospheric radii, while the mean radius of the maser 
spots remains at ~2.0 photospheric radii. 

Figure [12] shows a sketch of S Ori and its CSE near vi- 
sual minimum and post-maximum illustrating the findings of 
the present work. Our results on S Ori are largely consistent 
with the canonical scheme of the CSE of AGB stars discussed 
for instance in Reid & Menten ( I1997I I. However, our results in- 
dicate an inner dust shell boundary that is located closer to the 
optical photosphere compared to the Reid & Menten scheme at 
about two photospheric radii, and SiO maser shells that are co- 
located with the low-intensity molecular atmosphere and with 
AI2O3 dust grains near visual minimum. In the case of S Ori, 
the dust shell consists only of AI2O3 grains. For Mira stars that 
also show a dust shell of silicates, the inner boundary of silicate 
grains is still expected at greater distances of ^ 4 photospheric 
radii (Lorenz-Martins & Pompei a 120001 ). In addition to a static 
picture of a Mira star's CSE, our observations suggest the for- 
mation of a dust shell close to the surface near minimum phase 
and an expansion by ~ 20% of the inner dust-shell boundary 
toward post-maximum phase. 

Our modeling approach of adding the CLVs of a dust-free 
M model to that of a radiative transfer model of the pure dust 
shell includes two uncertainties. Firstly, the M model series is 
designed for the prototype Miras o Cet and R Leo and used 
here for the longer-period Mira S Ori. Secondly, we omit a pos- 
sible interaction of gas and dust outward of the inner dust shell 
boundary. Both effects might add some additional uncertainty 
to the absolute values of the obtained photospheric diameter 
values and dust shell characteristics. However, our data can be 
described well by our modeling approach and these additional 
uncertainties are not expected to be significantly larger than the 
diameter errors given in Table [T] and the step sizes of the dust 
shell parameters in Table H) Furthermore, the relative accuracy 
of the obtained parameters among the four epochs of obser- 
vations are clearly less affected by possible systematic uncer- 
tainties. Currently, the chosen approach appears to be the best 
available one. Similar studies of S Ori and other long-period 
Mira stars in the future would benefit from new dynamic model 
atmosphere series tailored to Mira stars with a longer period 
and a higher mass compared to o Cet and R Leo, as well as 
from self-consistent calculations of dynamic atmospheres and 
dust formation. 

Our mid-infrared interferometric data, in particular at the 
lower edge of MlDl's wavelength range, together with the use 
of the M model series provide constraints on the continuum 
photospheric angular diameter of S Ori. The resulting values 



are in good agreement with several earlier diameter measure- 
ments obtained in the near-infrared ^T-band by different instru- 
ments and at different stellar phases, which increases the confi- 
dence in our modeling approach. Additional confidence in the 
continuum photospheric diameters and the validity of the dy- 
namic model atmospheres can in future be reached by time 
series of high spectral resolution (R = 1000 - 10000) near- 
infrared interferometry with the VLTl/AMBER (Petrov et al. 
120031) instrument concurrent to VLTl/MlDl observations. 

Our results obtained at 4 epochs within one pulsation cy- 
cle suggest the formation of a dust layer consisting of AI2O3 
grains close to the stellar surface near stellar minimum and an 
expanded dust shell after stellar maximum. Our understanding 
of the connection between the dust formation frequency and 
the stellar pulsation will be further improved by using more 
epochs within one pulsation cycle. Dust-forming episodes may 
also have their own timescale that may not be directly related 
to the visual variability (e.g., Hofner & Dorfi 1997 ; Nowotny et 
al. 120051 ). a study of which requires observations extended over 
more than one pulsation cycle. Since the data of S Ori do not 
show any sign of silicate grains, similar observations of targets 
that belong to the intermediate class group of Lorenz-Martins 
& Pompeia (exhibiting AI2O3 and silicate grains) would be re- 
quired to study the detailed relationship among the AI2O3 shell, 
the silicate shell, and the photospheric pulsation. 

Both our MIDI data and the images of the SiO maser shell 
do not show any signature of a globally asymmetric distribution 
of dust or gas in the case of S Ori within our uncertainties. In 
general, asymmetric distributions of gas and dust are to be ex- 
pected for some Mira stars (e.g., Weigelt et al. il996i Boboltz 
& Diamond l2005l Tatebe et al. 120061) . A detection of asym- 
metric structures can be performed at the VLTI by sequentially 
using baselines of similar length but different orientation close 
in time/stellar phase to MIDI and by using the closure-phase 
instrument AMBER with baseline triangles comprised of base- 
lines with different orientation. 

The relative positions of star and SiO maser shells in our 
current study is unknown. Our comparison of the extension of 
the SiO maser shells to the continuum photospheric diameter is 
based on the assumption that the center of the maser spot dis- 
tribution coincides with the center of the star. This assumption 
is the most likely scenario for stars that do not show any sign of 
asymmetric dust and gas distributions, such as S Ori. In prin- 
ciple, the relative position of the maser shells could be related 
relative to the radio photosphere, which was detected by Reid 
& Menten ( 119971 ) for a sample of 6 Mira and semi -regular vari- 
ables. Additional measurements of H2O and OH masers rela- 
tive to the SiO masers could give additional information on the 
structure and kinematics of the more extended CSE. 
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Fig. 7. The (red) v — 2, J — 1-0 (42.8 GHz) and (green) v - I, J - 1-0 (43. 1 GHz) maser images overlaid onto pseudo-color 
representations of the infrared intensity. The continuum photosphere (in fact mostly hidden behind the molecular atmosphere in 
the A^-band) is enhanced to a light blue color. The darker blue shades represent our model intensity profile as in panels (e) of 
Figs-UHll and the green shades represent the location of the AI2O3 dust shell on top of the low-intensity extended wings of the 
molecular atmosphere. The true location of the star relative to the maser images is unknown. Here, we assume that the center of 
the star coincides with the center of the maser spot distribution. Synthesized beam sizes for epoch A are 0.59 x 0.40 mas and 
0.60x0.37 mas for the v = 1 and v = 2 transitions, respectively. For epoch B these sizes are 0.52 x 0.21 mas and 0.51 x0.20 mas, 
and for epoch C 0.46 x 0.16 mas and 0.44 x 0.16 mas, respectively. Epoch D is a MIDI-only epoch and SiO maser observations 
have not been obtained at this epoch. The size of each panel is 30 x 30 mas, corresponding to 14.4 x 14.4 AU. 
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Fig. 8. LOS velocity structure of the (top) v = 1,/ = 1-0 (43.1GHz) and (bottom) v = 2,7 = 1-0 (42.8 GHz) SiO maser 
emission toward S Ori. The epochs of observations are from left to right Epoch A (Oyis - 0.46), Epoch B (<l>vis = 0.56), and 
Epoch C (Ovis - 1.14). The top sub-panels of each epoch and transition show the spectra formed by plotting maser intensity 
versus velocity, color-coded in 1.7 km s"' velocity increments from redward to blueward. The solid lines in the top sub-panels 
represent the scalar-averaged cross-power spectra averaged over all of the VLBA antennas. The bottom (main) sub-panels plot 
the spatial and velocity distributions of the masers. The color of each point represents the coiTesponding velocity bin in the 
spectrum, and the size of each point is proportional to the logarithm of the flux density. Errors in the positions of the features are 
smaller than the data points. The dashed circles are based on the mean angular distances of the SiO masers from the centers of the 
distributions. The light blue filled circles in the centers illustrate for comparison the angular sizes of the continuum photosphere 
derived in Sect. 13.41 The true location of the star relative to the masers is unknown and here assumed to coincide with the center 
of the maser spot distribution. We use also panels (e) of Figs.llHl] as well as Fig.|7] for a comparison of the SiO maser ring radii 
to the extension of the molecular atmosphere and the inner dust sheU boundary. 
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S Ori, near visual minimum 
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Fig. 12. Sketch of the radial structure of S Ori's CSE at (left) near-minimum and (right) post-maximum visual phase as derived in 
this work. Shown are the locations of the continuum photosphere (dark gray), the at A^-band optically thick molecular atmosphere 
(medium dark gray), the at A^-band optically thin molecular atmosphere (light gray), the AI2O3 dust shell (dashed arcs), and the 
42.8 GHz and 43.1 GHz maser spots (circles/triangles). The numbers below and beside the panels are the mean values of (left) 
epochs A & B and (right) epochs C & D from Table |2l 



